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Neutron resonance spin-echo spectroscopy was used to monitor the temperature evolution of the 
linewidths of transverse acoustic phonons in lead across the superconducting transition temperature, 
Tc, over an extended range of the Brillouin zone. For phonons with energies below the supercon- 
ducting energy gap, a linewidth reduction of maximum amplitude ~ 6^eV was observed below Tc- 
The electron-phonon contribution to the phonon lifetime extracted from these data is in satisfactory 
overall agreement with ab-imtio lattice-dynamical calculations, but significant deviations are found. 
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The electron-phonon interaction is the major limiting 
factor for electronic transport phenomena in metals at 
elevated temperatures, and it is responsible for most in- 
stances of superconductivity. For metallic and supercon- 
ducting elements and binary compounds such as MgB2 
0, modern ab initio calculations yield accurate predic- 
tions for the electron-phonon coupling parameters of ev- 
ery phonon over the entire Brillouin zone (BZ). An ac- 
tive area of research is aimed at a realistic description of 
electron-lattice interactions in complex compounds with 
strong electronic correlations. Yet only rudimentary ex- 
perimental tests of these calculations have thus far proven 
possible. For instance, electronic transport or tunnelling 
experiments probe weighted averages of the electron- 
phonon interaction over the entire phonon spectrum. Op- 
tical spectroscopy is capable of probing the lifetimes of 
individual phonons limited by scattering from electrons, 
but kinematics constrains these experiments to a single 
point in momentum space. Here we use a new neutron 
spectroscopy method |M ^ IS with an energy resolution 
in the /zeV range (that is, about two orders of magnitude 
better than that of standard neutron spectroscopy) to 
determine the electron-phonon lifetime of an individual 
acoustic phonon in lead over an extended range of the BZ. 
The results are compared to ab initio lattice dynamical 
calculations also reported here. Our experiment consti- 
tutes the first detailed test of modern calculations of the 
electron-phonon interaction and opens up a new avenue 
for a quantitative understanding of electron-phonon in- 
teractions and superconductivity in solids. 

Since the 1950's, triple-axis spectrometry (TAS) with 
neutrons (and recently, x-rays) has been the method 
of choice to experimentally determine energy- and 
momentum- resolved phonon spectra of solids. Briefly, 
the energies and momenta of the incoming and scat- 
tered neutrons are selected by crystal monochromators, 
and their difference yields the phonon dispersion rela- 
tion. The monochromaticity of the beam is thus coupled 



to the beam divergence, which can only be restricted at 
the expense of beam intensity. This implies that energy 
resolutions significantly better than 10% are impracti- 
cal to achieve under almost all circumstances. Since the 
electron-phonon interaction leads to a phonon linewidth 
broadening of typically less than 1%, TAS allows phonon 
lifetime measurements only in exceptional cases. While 
the phonon dispersions in lead were determined by TAS 
early on , attempts to resolve the electron-phonon life- 
times have been unsuccessful 0, 1^ • TAS is therefore not 
generally applicable as a probe of the electron-phonon 
interaction in metals. 

A neutron spin-echo method to improve the energy 
resolution of TAS by several orders of magnitude with- 
out loss of intensity was proposed some time ago 21 ■ 
The method has been described in detail elsewhere ^. 
Briefly, it is implemented by inserting spin polarizers and 
tunable Larmor precession coils into the incident and 
scattered beams of the triple-axis spectrometer. If the 
phonon to be measured has infinite lifetime, and if the 
precession fields are properly matched (spin-echo condi- 
tion), the full polarization of the incident beam is re- 
covered at the detector. If the phonon lifetime is fi- 
nite, the spin-echo amplitude is reduced, and the phonon 
linewidth can be determined by systematically varying 
the precession fields. Early versions of this technique 
have enabled the determination of the lifetimes of ro- 
tons in superfluid helium and optical phonons near 
the BZ boundary in germanium However, a major 
technical obstacle has thus far prevented more general 
applications. The dispersive nature of most collective 
excitations in solids leads to a degradation of the neu- 
tron spin polarization that masks the effect of the finite 
excitation lifetime. This effect can be compensated to 
linear order by tilting the precession coils away from the 
neutron beam direction by an angle proportional to the 
slope of the dispersion relation . Since the required tilt 
angles are of order 10-50°, however, this cannot be ac- 
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complished by the long solenoids used in the early work. 

In the "neutron resonance spin-echo" (NRSE) tech- 
nique H, ^ 1^, the solenoids are replaced by a pair of 
compact radio-frequency (r/) coils surrounding a field- 
free region. This allows larger tilt angles sufficient to 
match the dispersion relations of most collective excita- 
tions in solids. Using a prototype NRSE-TAS setup, the 
feasibility of phonon lifetime measurements in lead was 
demonstrated at temperatures exceeding 50 K, where the 
lifetime is limited by phonon anharmonicity |l2| . How- 
ever, the neutron flux was insufficient to obtain high- 
quality data at lower temperatures, where the electron- 
phonon interaction is dominant. 

Here we report the results of experiments carried out at 
the dedicated, high-flux NRSE-TAS spectrometer TRISP 
at the FRM-H neutron source in Garching, Germany 
Starting from the source, in-beam components 
include a spin-polarizing guide; a velocity selector to 
cut out higher-order contamination of the incident beam; 
a horizontally and vertically focusing pyrolytic-graphite 
(PG) monochromator and a horizontally focusing PG an- 
alyzer, each set for the (002) reflection; and a supermirror 
spin-polarizer in front of the detector. The measurements 
closely followed the protocol established in Ref. ^3- The 
reader is advised to consult this reference for a detailed 
description of the experimental setup and data analysis. 
The samples were two cylindrical Pb single crystals with 
cylinder axes (100) and (110) and mosaicities 4.2' and 
5.7', respectively. Both crystals had a diameter of 30 
mm and a length of 50 mm. The first crystal was used 
to measure the Ti-phonon along (110) and the T-phonon 
along (100) in the scattering plane spanned by the (001) 
and (010) directions of the reciprocal lattice. The re- 
maining measurements were carried out with the second 
sample. The samples were loaded into a closed-cycle He 
cryostat. Depending on phonon energy, the incident (fi- 
nal) neutron wave vector was fixed at 2.5lA~^ (1.7A~^). 
The triple-axis spectrometer was aligned in a defocusing 
configuration in order to minimize the segment of the 
dispersion surface sampled by the TAS resolution ellip- 
soid. The spectrometer angles were tuned to maximize 
the phonon intensity and kept constant during the spin- 
echo scans. The NRSE coils between monochromator and 
sample and between sample and analyzer were tilted to 
match the slope of the acoustic phonon dispersions. The 
spin-echo time, r, was changed by varying the rf frequen- 
cies in the NRSE coils. In order to minimize systematic 
errors, the beam polarization at the detector, P, was de- 
termined by translating the coil immediately in front of 
the detector and extracting the amplitude of the result- 
ing sinusoidal intensity modulation from a least-squares 
fit [13. 

The strategy of the experiment was to single out 
the electron-phonon contribution to the phonon life- 
time by monitoring the evolution of the phonon line- 
shape through the superconducting transition tempera- 
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FIG. 1: (Color online) Neutron resonance spin echo profiles 
for the Ti-phonon at (0.3, 0.3, 0) at selected temperatures 
above and below the superconducting transition temperature. 
The inset shows a triple-axis scan in focusing mode for the 
same phonon. From these data, F was extracted as 25.5 ±0.6, 
31.4 ± 0.8, and 30.9 ± 0.7 fieV for T = 3.7, 5.5, and 8.0 K, 
respectively. 



ture, Tc ~ 7.2 K. As demonstrated previously by TAS for 
phonons with anomalously large linewidths in NbsSn [l5| 
and Nb [l^ , this contribution vanishes if the phonon en- 
ergy is below the superconducting energy gap, 2A. Other 
intrinsic contributions to the phonon linewidth, such as 
isotope disorder and the anharmonicity of the lattice po- 
tential, do not exhibit anomalies at Tc. 

Because of its face-centered cubic lattice structure, the 
phonon spectrum of lead consists entirely of acoustic 
branches. Most of the data were taken on the lowest-lying 
transverse mode, Ti, in the (^^0) direction, where all 
phonons are nondegenerate. The atomic displacements 
of this mode are along (110). The mode energy is lower 
than 2A over much of the BZ, so that a large range of 
phonon momenta can be covered in this way. More lim- 
ited data sets were also acquired for the T2 mode (with 
displacements in the (001) direction) along (^^0), and for 
the transverse phonon along (^^^). 

Fig. 1 shows typical NRSE data for the Ti mode 
along (C^O). The dependence of the beam polarization on 
the spin-echo time, P(r), is proportional to the Fourier 
transform of the phonon lineshape 0. As a Lorentzian 
lineshape thus yields an exponential NRSE profile, the 
phonon linewidth can be obtained from the slope, F, of 
(P(t) on a semilogarithmic plot (lines in Fig. 1). A 
drop of the slope (corresponding to an increase of the 
phonon lifetime) below Tc is clearly apparent. In addi- 
tion to the intrinsic phonon linewidth, F is also affected 
by extrinsic factors such as the mosaicity of the crystal 
and the curvature of the dispersion relation, as discussed 
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FIG. 2: (Color online) a-e) Temperature dependence of the 
observed NRSE relaxation rate F for Ti-phonons alone (^^0). 
f) Superconducting energy gap from tunnelling data p^ . g,h) 
F versus temperature for r2-plionons along (^^0), and i) for a 
transverse phonon along (5^0- The dashed lines in panels a-e, 
g, i indicate the temperatures at which the superconducting 
gap opens, according to the data of panel f. 



in detail in Refs. However, these temperature- 

independent factors have a negligible influence on the 
linewidth anomaly at Tc, which is of interest here. A 
detailed discussion of resolution effects on the present 
data will thus be given elsewhere the following 

figures presenting the experimentally observed F versus 
temperature, the extraneous contributions give rise to 
zero-temperature offsets that generally decrease with in- 
creasing wave vector, due to the decreasing curvature of 
the dispersion surface sampled by the TAS resolution vol- 
ume. 

Fig. 2a-e shows the temperature dependence of F for 
selected Ti-phonons along (•C'^O). A decrease of F at Tc is 
evident in all of the traces shown for < 0.32 {E < 2.46 
meV). The amplitude of the effect increases with increas- 
ing ^ and reaches ~ 6/ieV for ^ = 0.32. For ^ = 0.4 
{E — 3.05 meV), the effect is no longer present within 
the experimental error. This observation is in good ac- 
cord with the magnitude of the superconducting energy 



gap, 2A 2.6 meV, extracted from tunnelling data [l^ . 
For £^ = 0.43, there is some indication of an increase in 
phonon linewidth at low temperatures. This is expected 
because of the pileup of electronic states above the energy 
gap. This effect is also responsible for the nonmonotonic 
temperature dependence of F for ^ = 0.32 (Fig. 2c). As 
2A approaches the phonon energy E from below with de- 
creasing temperature, F first increases gradually because 
of the enhanced density of states above the gap, and sub- 
sequently plummets sharply below Tc as 2A > E. It has 
been shown [lH.ll9|| that the jump in phonon linewidth at 
Tc, expressed as a fraction of the normal-state linewidth, 
is equal to the jump in the ultrasound attenuation co- 
efficient at Tc- This in turn can be calculated in BCS 
theory [l^ as (7r/2)[l - 2/(1 + cxp(£;/fcBT))] = 1.44 for 
E 2.46 meV and T = 4.6 K (the temperature where 
2A — E), in good agreement with the value 1.48 ± 0.1 
determined from Fig. 2c. 

For the more steeply dispersing T2 phonon along (i^i^O), 
the superconductivity-induced phonon linewidth renor- 
malization is much smaller than that of the Ti phonon, 
and only upper bounds could be established (Fig. 2g). 
This is expected on general grounds, because the steep 
dispersion of this mode implies that the linewidth effects 
are restricted to wave vectors close to the BZ center. In 
the (^^^) direction, where the two transverse phonons 
are degenerate, we observed a superconductivity-induced 
linewidth reduction of magnitude similar to that of the 
Ti-phonon along (^^0) (Fig. 2i). 

In order to compare the experimental results to predic- 
tions of ab-initio calculations, the phonon dispersions and 
linewidths were calculated in the framework of density 
functional [13, 0| perturbation theory [i^ in the local- 
density approximation. We employed ultra-soft pseu- 
dopotentials ,23] , with scalar-relativistic corrections and 
a plane- wave basis set |2^, with a cut-off energy of 40 
Ryd. The k-space integration was approximated with a 
(12)^^ Monkhorst-Pack grid [2^ using a smearing param- 
eter [2^ of 0.06 Ryd for the self-consistent cycles and 
phonon calculations, and with a much denser (36)^^ mesh 
and a Gaussian smearing of 0.04 Ryd for the calculation 
of the linewidths. With these parameters and the numer- 
ically optimized lattice parameter of 4.95A, the phonon 
frequencies were converged to 0.1 meV and the linewidths 
to 0.025 ^eV. (If the lattice parameter 4.92A determined 
experimentally at low temperature is used instead, the re- 
sults are not affected outside these numerical confidence 
limits.) The calculated phonon frequencies agree with 
the experimentally determined phonon frequencies 0, Q 
within the published error bars. For selected k-points, 
we recalculated the phonon energies and linewidths using 
Savrasov's full potential LMTO program ^ and found 
very good agreement. 

Fig. 4 provides a synopsis of the experimentally 
determined and numerically calculated electron-phonon 
linewidths along (^^0). Although the overall agreement 
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FIG. 3: (Color online) Experimentally determined and nu- 
merically calculated electron-phonon linewidth of transverse 
phonons along (^^0)- The experimental points were obtained 
by subtracting the linewidths at lOK (see Fig. 2) from those 
at the lowest temperature (0.5 K). 



is satisfactory, a significant, systematic deviation outside 
of the experimental error bar is apparent around ^ = 0.3. 
A similar deviation was found for the more limited data 
set along {^£,£,). A possible source of this discrepancy is 
the spin-orbit coupling, which was not considered in the 
calculation. The electronic density of states at the Fermi 
level increases by about 10% when spin-orbit coupling 
is incorporated in the electronic structure calculation. A 
corresponding increase in the electron-phonon interaction 
is to be expected. Indeed, Pb is a well-known example 
in which the agreement between the experimentally de- 
termined and theoretically calculated phonon frequencies 
and total electron-phonon coupling parameter A is not as 
good as in the case of other elemental superconductors, 
such as Al or Nb. In our calculation, A was computed 
as 1.27, whereas A = 1.55 was determined by tunneling 
spectroscopy 0. Furthermore, the values of A calcu- 
lated from first principles range from A = 1.2 (Ref. l27j) 
to A = 1.68 (Ref. J^). The experimental results presented 
here therefore call for further development in the numer- 
ical computation of electron-phonon interactions. 

In conclusion, we have used neutron resonance spin- 
echo spectroscopy 0, 0, to determine the electron- 
phonon contribution to the lifetimes of selected phonons 
in lead. The results have enabled the first energy- and 
momentum-resolved test of ab initio lattice-dynamical 
calculations. The energy resolution of the experiment 
was about two orders of magnitude better than that 
achievable by triple-axis spectroscopy (for a typical 
triple-axis scan, see the inset in Fig. 1). The method 
can be applied in a straightforward manner to other el- 



ements and compounds that exhibit superconductivity 
at low temperatures. For non-superconducting metals, a 
rigorous treatment of the resolution function is required 
to extract the electron-phonon lifetime. The pertinent 
formalism has already been established [l^ . Finally, 
the technique has recently been successfully applied to 
magnetic excitations ^ and opens many new perspec- 
tives in this field as well. 
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preparing the Pb crystals. 
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